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Abstract 

Heavy nuclei at collider energies are a source of high energy Weizsacker-Wihiams photons. 
This photon flux may be utihzed to study high energy photon-nucleus interactions. Here we 
concentrate on the coherent diffractive production of heavy vector mesons on nuclear targets 
and show how it probes the unintegrated glue of the nucleus in the saturation domain. We 
present predictions for rapidity distributions of exclusive coherent J/^ and T mesons which 
can be measured by the ALICE experiment at the LHC. 
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I. INTRODUCTION 



The exclusive photo- and electroproduction of vector mesons at high energies has re- 
cently been thoroughly studied at the electron-proton collider HERA (for a review, see 
In a Regge picture, this process is driven by a t-channel Pomeron exchange. The 
HERA data, which range from low to high photon virtuality as well as from light to 
heavy vector mesons, have given an intriguing insight into the Pomeron physics from 
soft to hard processes. Overall a consistent phenomenology emerges, in which the QCD- 
Pomeron exchange is modelled by the gluon-ladder exchange and quantified by the unin- 
tegrated gluon distribution of the target proton or the color dipole-proton cross section. 

In recent years, the very-small x behaviour of the gluon structure function has been 
of great interest in the context of unitarity effects and the saturation phenomena 

The multiple scattering/absorption effects associated with the gluon saturation 
physics are naturally enhanced on a large nuclear target. From this point of view 
the coherent diffractive production of vector mesons on heavy nuclei is very interest- 
ing. Clearly, the best option to study the small-x nuclear glue would be a dedicated 
electron- ion collider jsl, where for example one could measure nuclear structure func- 
tions at perturbatively large Q^. Lacking such a facility, there appears to be no easy 
experimental access to the nuclear gluon distribution. 

Here, the exclusive photoproduction of heavy vector mesons J/\E', T opens up a new 
possibility to experimentally study the small-x nuclear glue. It is the large mass of the 
heavy quark which provides the hard scale and assures that we may use perturbation 
theory even in the photoproduction limit. 

Presently, nucleus- nucleus collisions at ^/snn = 2.76 TeV studied at the LHC offer 
an access to high-energy photonuclear reactions [4]. Here, one of the nuclei will serve 
as a source of Weizsacker- Williams photons, while the other one plays the role of a 
target. Coherent diffractive processes of interest in this work leave the target intact and 



deflected only by a very sma^ 
on nuclei is an old subject 



I angle. Diffractive photoproduction of light vector mesons 
si, and with some success can be addressed using vector- 
meson dominance, and within a limited range of energies, a hadronic Glauber coupled 
channel type of model. When a hard scale is present - such as a large photon virtuality. 
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or a large quark mass - it proves more efficient to use perturbative QCD degrees of 
freedom. In particular the color-dipole formulation to which we will now turn allows to 
account for nuclear effects in a rather straightforward manner. 



II. AMPLITUDE AND CROSS SECTION FOR ^ VA REACTION 
A. Vector meson production in the color-dipole picture 
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FIG. 1: Sample diagram for the — t- VA amplitude. At x ~ the relevant degree of freedom 
is the QQ dipole. All multipomeron exchange diagrams are summed up by the Glauber- Gribov 
multiple scattering theory for color dipoles. 



It is useful to start from the color-dipole formulation {g] of vector meson photoproduc- 
tion. This formulation is essentially equivalent to the A;^-factorization approach which 
we previously used for the free proton target [Tj], and to which we turn later in this work. 
The multiple scattering corrections relevant for nuclear targets are however more easily 
derived in impact parameter space. 

The forward amplitude for vector-meson photoproduction on the nuclear target at 
center-of mass energy W (x = rriy/W'^), takes the form 



Ai^A^Vp;W) = t{V\aAix,r)\^) = 2t J d^b{V\TAix,b,r)\^) 
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2iy" d'^b J dz J d'^r tpy{z,r)ip^{z,r)TA{x,b,r) . 



Here iljv{z,r),tlj^{x,r) are the light-cone wave-functions of the vector meson and 
the photon, respectively. For simplicity we suppress the summation over (anti-)quark 
helicities. The explicit form of the wave function are not important for the argument in 
this section, we will later give all explicit formulas in momentum space. 

The amplitude can be easily generalized to finite transverse momentum transfer A: 

A{-fA Vp;W,A) = 2i j d^b expH^A] {V\V A{x,b,r)\^) . (2) 
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The differential cross section is then given by (t = —A 



dt dA' 167r 



Following 



at small x {x < xa = O.IA^^^^), the multiple scatterings of the color- 
dipole can be summed up by a Glauber-series for color dipoles. The color-dipole-nucleus 
amplitude in impact parameter space for x ^ xa can then be given in terms of the 
color-dipole-proton cross section: 

Ta{xa, b, r) = 1 - exp ( - ^a{xA, r)TA{b)) . (4) 

In practice that means xa ~ 0.01. For smaller x one must take into account higher, 
qqg-Fock states as we discuss below. 

In order to quantify the size of nuclear (multiple scattering) effects, one often compares 
to the impulse approximation. The latter works well if multiple scatterings are weak, i.e. 
if the nuclear opacity a{xA, r)TA{b)/2 is small. In impulse approximation we assume that 
only one of the nucleons in the nucleus participates in the interaction, and all others are 
spectators. Expanding the Glauber-exponential to the first order we obtain (hereafter 
lA stands for impulse approximation): 

Aia{iA VA; W,A)=i {V\a{x,r)\-f) J d'^bexp[-tbA]TA{b) . (5) 

The total cross section in impulse approximation would then be 

atot;iA(7A ^ VA; W) = 4 ^^^^^ / d^TKb) . (6) 
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Let us introduce the ratio of the full nuclear cross section to the impulse approximation 
result: 



fTtot,iA(7^ ^ VA] W) 
Then we can also express the total photoproduction cross section on the nucleus as: 

daijp — > Vp) 



(7) 



atot(7^ ^ VA; W) = i?eoh(W^) ■ 47r- 



dt 



t=o 



Jd'bTlib). (8) 



Here 



da{^p Vp) 



dt 



t=o 



By ■ (TtotilP Vp) , 



(9) 



where By is the diffraction slope, can be taken from experimental data 



Hii|. 



le integral over the nuclear optical density squared behaves parametrically as (see 



e.g. 



12|): 



jd'bTl{h) = CA 



4^' 



(10) 



where Rch is the nuclear charge radius, and Ca is a number of order unity which depends 
on the shape of Ta (b) . In the numerical calculations, we use a realistic nuclear density. 



as parametrized in 



13|. 



B. Momentum space formulation of vector- meson production on nuclei 



We can bring the photoproduction amplitude for the nuclear target into the similar 
/c^-factorization form as the result for the free proton. The only difference is that now 
the unintegrated gluon distribution of the proton, will be replaced by the appropriatly 
defined unintegrated glue of the nucleus, explicitly constructed in the treatment of the 
diffractive tcA — )■ qqA process in |l2 |. 

Recall now, that for the free-nucleon target, dipole cross section and unintegrated 
gluon distribution are related by p^] 

a{x, r) = (To / d'^K, 
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where we pulled out ctq = cq^x), so that / is normalized to unity: 



/(x,k) = — — — 



ao N, 91og(K2) 



CTo(x) 



47r 1 dGN{x,H? 



k4 91og(K2) 



(12) 



We can analogously introduce the impact-parameter dependent unintegrated gluon dis- 
tribution of the nucleus Mb, ) through the relation 



Ta{x, b,r) = / (fn 



1 - e 



iH,r 



[b,X, Av) . 



(13) 



(14) 



In terms of the nuclear glue Ga, the function (f){b, x, n) fulfills 

nO,X,H.)- d log K^d^b • 

n the forward scattering limit, i.e. for A = the photoproduction amplitude given in 
l| can be brought in the form (see 7]): 

1 dz 



3mr(iy, = 0) = |^2C.V47rae™ 



47r2 



z(l — z) Jo 



Txdk ^viz, k 



(15) 



K"^ 



asiq^J^Aix, k') iAo{z, k') Wo{k\ k') + A^{z, k') Wi{k\ k') 



where 



Aoiz,P) = ml + 



k'^m, 



Q 



M + 2mQ 



Ai{z,e) = z^ + {1- zf -{2z-l) 



2 ruQ 



M + 2mQi P + m| 



,2 ' 



Wo{k',K') 



1 



k"^ + m% 



Wiik^,K^'' 



k"^ + ml 



^Jik"^ — ttlq — K^Y + ArnQk"^ 
k!^ — — 



Q 



2P 



^{k"^ — rriQ — K^Y + Arnqk"^ 



(16) 
(17) 

(18) 
(19) 



Here mq is the heavy quark mass, and 



„ k"^ + m% , ^ 

= — ^ 20 

z{l-z) 

is the invariant mass squared of the QQ system in the final state. The strong coupling 
enters at the hard scale = max(K^, k"^ + rn^). The light-cone wave function tpy of the 
vector meson is parametrised exactly as in [l|, l7|, and the Gaussian form, which proved 
to lead to good agreement with experiment, is adopted. 
All nuclear effects are accounted for by the substitution 

asdGA{b,x,Hi'^) 



^^J^a{x,^') = Jd'b 



d log K'^d'^b 2n 



d b(f){b, X, At) 



(21) 



In this way we introduce also the impact-parameter dependent amphtude by the relation 

Tir^A ^VA) = J (fb QmT{-fA VA; b) . (22) 



Its normalization is such that 

da{-fA -> VA; W) 1 ^niTi-fA -> VA; b) 2 



d?b 4 



(23) 



C. The nuclear unintegrated glue in the Glauber regime 

Here we briefly recapitulate how to calculate the nuclear unintegrated gluon distribu- 
tion from the proton unintegrated glue. 

The starting point is the definition of the nuclear unintegrated glue (fT3|) . We are 
interested in the regime of a; ~ x^, where the Glauber representation of the color-dipole 
scattering amplitude r(x, b, r) is valid: 

T{xA, 6, r) = 1 - exp[ia(xA, r)TA{b)] . (24) 

Introducing the shorthand notation 

v{b) = ]^as ao(x^) TA{b) , (25) 

and the the multiple convolutions 

/(^■)(x, K.) = j d'K, . . . d'K,6^'\K - Y: >^i) fix, K,)... fix, K,) , /(°)(a.) ^ 6^'\k) , (26) 

i 

we obtain the expansion of the Glauber exponential 

exp[-ior(xA,r)T^(6)] = ^ [ d^n. exp[iKr]wkib) f^'\xA, i^) (27) 
^ k>o •' 

with the Poisson-weights 

«;fc(b) =exp[-z/(b)]-^. (28) 

This gives us the expression of the nuclear unintegrated gluon distribution as an expan- 
sion over multiple convolutions of the free nucleon glue: 

<j,ib,x,K) = Y.w,ib)f^^\x,K). (29) 
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here, the j-th term of the expansion is the contribution of the interaction of j nucleons 
with the color dipole. Here f^^\x,K) is the collective unintegrated glue of j nucleons, 
and the weight factor Wj gives us the probability the j nucleons overlapping at impact 
parameter b take part in the interaction. 

It is important to realize that the nuclear unintegrated glue includes the multiple 
scattering corrections, it is not simply a two-gluon exchange in the crossed channel, as in 
the free nucleon glue. To some degree however, these multiple gluon exchanges behave 
like a two-gluon exchange: the diffractive amplitude has the same form as the two-gluon 
exchange amplitude on the free nucleon target. 

D. Small-x-evolution: contribution of the QQ^-Fock state 



9 




FIG. 2: A sample diagram which contains the multiple scattering of a QQg-Fock state. This 
is the first step in the nonlinear evolution of the nuclear unintegrated glue. 

If we increase the jA center-of-mass energy successively higher and higher 
QQg, QQgg, ■ ■ ■ Fock-states become important. On the free nucleon target, the effect of 
higher Fock-states, which contain gluons strongly ordered in rapidity can be resummed, 
and leads to the small-x evolution of the then x-dependent dipole cross section [l^. The 
multiple scatterings of the QQg-Fock state off a heavy nucleus are given by the first 
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iteration of the nonlinear Balitsky-Kovchegov 15| evolution equation 

^J'^^nf} = ^BFKL ® m ^, P) + ^, P) ■ 

alog^l/x) 

The nuclear glue, which includes rescattering corrections of QQ as well as QQg Fock 
states is then given by 

(f){b,x,p) = (f){b,x A, p) + log [ — ) ■ . 30 

In a similar manner, to obtain the ratio i?coh of eq.([7]), the impulse approximation am- 
plitude is calculated from 

<PMb,XA,p) = TAb) . - -^^-^ , (31) 

subject to a similar iteration ( |30|) as the full glue, but with the nonlinear piece omittec . 



For the explicit momentum-space form of the infrared-regularized BK-equation, see |16 |. 
A similar strategy of including the QQ and QQg-Yook states has been followed for 
the nuclear structure function and inclusive coherent diffraction in Ref. [l3]. There a 
good agreement with available data on nuclear shadowing has been obtained. In the 
numerical calculations, we use an unintegrated gluon distribution of the proton that has 



been obtained from an analysis of HERA structure function data in 



18|- 



Let us briefly discuss, how our approach differs from others available in the literature. 



The first estimates of Klein and Nystrand 



19| are based on extracting an effective J j'^- 



nucleon cross section from photoproduction data using vector-meson dominance ideas. 
They then go on to use this cross section to evaluate the classical survival probability of 
mesons passing through a slab of nuclear matter. Goncalves and Machado [20] adopt the 
color dipole approach and give a proper quantum mechanical treatment of the multiple 
scattering effects. Their approach differs from ours in that they absorb all saturation 
effects into the dipole-nucleon cross section which is then eikonalized. This is strictly 
speaking inconsistent with the nonlinear evolution of the dipole-nucleus cross section, 
and neglects the fact that multiple scatterings off different nucleons are enhanced by 
the nuclear size. Of course it may well be viable phenomenologically in a limited range 
of energies. Finally, Rebyakova et al. use a relation of the diffractive amplitude 
to the integrated gluon distribution of the target, which holds, with some reservations, 
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for heavy quarks. Such an approximation can be obtained from the /c^-factorization 
formahsm used in this work to the leading logarithm in the hard scale (see |l| and 
references therein). There appears to be a hidden assumption that all saturation effects 
are summed up in a boundary condition of the integrated, coUinear, nuclear glue. 

III. RESULTS AND CONCLUSIONS 

A. Photoproduction on nuclei: 'yA — )• VA 



^ 1 




0.6 - 




'—-—J ■ — ■ — ■ — I 

10 10^ 10^ 

W (GeV) 

FIG. 3: Ratio of the nuclear coherent cross section for J/^' and T production to the Impulse 
Approximation cross section. 

In Figj3]we show the ratio Rcoh{W) for the lead nucleus (^°^Pb) for J/\E' (red line) 
and T (blue line) meson production. The deviation of Rcoh from unity is a measure of 
the strength of nuclear rescattering/absorption effects. We see that the nuclear effects 
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are stronger for the J/\I/ than for the T meson. This is indeed to be expected, as we can 
most easily see in the dipole picture: the photoproduction amphtude probes the dipole 

n 

cross section at the scanning radius oc l/rriQ [6], and the smaller dipoles relevant for 
T-production will experience smaller rescattering effects. 

In FigJH we present the impact parameter distribution of vector meson photoproduc- 
tion. We show results for J/\l' and T mesons, for two different energies (W = 200, W = 
2760 GeV) and lead nuclei. 
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FIG. 4: Impact parameter distributions for J/\I' and T mesons for the ^"^^Pb target. 



B. Ultraperipheral nucleus-nucleus collisions: AA AAV 

In the AA — )■ AAV processes the heavy nuclei play two different roles. One of the 
nuclei is a target and the next is a source of high energy Weizacker-Wiliams photons. 
In Fig.(|5]) we show Feynman diagrams for the relevant production mechanism, the Born 
diagram in the left panel and a diagram including the absorptive correction in the right 
panel. It should be noted that one should also add the amplitude in which the photon 
is emitted from the lower line. The interference of the photon emission from the upper 
or lower lines in fact causes a peculiar azimuthal correlation between the outgoing nuclei 
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22| . After integration over azimuthal angles, and at the Born level, the interference 
drops out and we can add the squares of both amplitudes. If absorptive corrections are 



included, a small interference contribution remains even after azimuthal averaging |22|. 
Below we will neglect the interference effect and evaluate the nucleus-level cross section 
from the absorption corrected equivalent photon approximation: 

aiAiA^ ^ A1A2V; s) = J du 2^ aijA^ VA^] 2oj^s) + {I ^ 2) . (32) 

To obtain the effective photon flux dN^^, one starts from the electric field strength 
associated with the moving nucleus (see for example Ref. |4j for a review and references) 



/d^q 
j^—:^ exp[~ibq] 



(33) 



q^ + ixP- 1^'^ 

Here FemiQ"^) is the charge form factor of the nucleus, u is the photon energy, 7 the 
relativistic Lorentz-boost of the beam and b is the impact parameter. Having the electric 




FIG. 5: Mechanism of exclusive vector-meson production in AA collisions. Left panel: Born 
diagram; Right panel: the elastic rescattering in the initial state accounts for the unitarity 
effect of inelastic channels. 



field strength as a function of photon energy and impact parameter we can calculate the 
photon flux corresponding to the Born diagram in the left panel of Fig. (5) 



dN{u, b) 



dud'^b ^, ,,2 
E{uj, b) 



(34) 
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Finally, the effective photon ffux with absorptive corrections included is: 

dN^^f{uj) = J d^bSli{b)dN{u,b) . 

The absorptive correction S'^i{b) is shown schematically by the extra oval in the right 
panel of Fig. O These absorptive correction can be calculated by applying the following 
simple formula (see e.g. [23|): 



S^iib) = exp ( - aMNTA.AM) ~ ^(l&l " (^i + ^2)) 



(35) 



where Ri and R2 are the radii of the colliding nuclei. We remove those configurations 
in the impact parameter space, when the nuclei overlap, which at high energy means 
automatically their break up. Therefore absorptive corrections have a meaning of the 
gap survival probability. In Fig. E] we present the total cross section for J/\l' and T 
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FIG. 6: Total cross section for Pb Pb Pb J/* Pb and Pb Pb ^ Pb T Pb as a function of 
nucleon-nucleon energy. 

production in nucleus-nucleus collisions as a function of energy for the ^osp^ target. The 
cross section for J/\I' production increases by two-orders of magnitude when going from 
RHIC to LHC energy. For the T meson the corresponding increase of the cross section 
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FIG. 7: Rapidity distribution of J/^ for symmetric collisions of lead nuclei for Wmn = 200 
GeV (upper) and Wat at = 2760 GeV (lower). Individual contributions are shown separately. 

is substantially bigger. In Figs. [7] -[8] we show the differential cross section in rapidity 
for the exclusive coherent productionof J/\l/ and T mesons, in lead- lead collisions. The 
shape of the distributions strongly depends on the collision energy. The distributions in 
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FIG. 8: Rapidity distribution of T for symmetric collisions of lead nuclei for Wnn = 200 GeV 
(upper) and Wnn = 2760 GeV (lower). Individual contributions are shown separately. 

rapidity for AA collision are much narrower than similar distributions in proton-proton 
collisions [3, Q]. This has its origin in the large nuclear size: the charge form factor 
F{q^) of the nucleus is much sharper compared to the proton's Dirac form factor Fi{q^). 
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The spectrum of Weizsacker- Williams photons in a proton is considerably harder than 
in the nucleus, see e.g. {4]. 

In summary, we presented predictions for the exclusive coherent diffractive production 
of J/^ and T-mesons in collisions of heavy nuclei at LHC energies. Our framework takes 
into account not only the Glauber-type rescattering of color dipoles in the nuclear matter, 
but also the gluon-fusion/shadowing corrections associated with the rescattering of the 
Q(Q(7-Fock-state. The predicted rapidity distributions of mesons may be tested by the 
ALICE experiment at LHC 
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